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Abstract 

We apply an algebraic method for studying the stability with respect to a set of conserved 
quantities for the problem of torque-free gyrostat. If the conditions of this algebraic method are 
not fulfilled then the Lyapunov stability cannot be decided using the specified set of conserved 
quantities. 

! MSC 2010: 34D20, 37B25, 70E50, 70H14 

"■j"^' Keywords: stability, rigid body, gyrostat 

1 Introduction 

A very important problem in the theory of the differential equations is the problem of the stability A 
very useful tool for determining stability of an equilibrium point is Lyapunov's direct method connected 
^ ' with the Lyapunov functions. A natural candidate to be a Lyapunov function is a conserved quantity 

In a lot of examples coming from mathematical physics, we identify a set {Fx, Fk} of conserved 
quantities. In many situations they are not positive definite functions in the equilibrium points of 
interest. In this situation, a first step to decide if the equilibrium point is stable is to search a Lyapunov 
function of the form $>(Fx, ...,Fk), where $ : M fc — > R is a smooth function. A function $(Fi, i^) 
is a Lyapunov function if and only if it is a positive definite function. In Stability Theory are known 
some important methods to construct positive definite functions using conserved quantities. We remind 
the so-called "Chetaev's method" presented in [TT] and some methods which appeared in the context of 
Hamilton-Poisson systems. In 1965 have been introduced, see [2] the Arnold's method. At the beginning 
of eighties was developed the Energy-Casimir method (see [S], [B]) and in 1998 the paper [7] present 
the Ortega-Ratiu method. In [3] is proved the equivalence of the Arnold's method, the Energy-Casimir 
method and the Ortega-Ratiu method. 

If there exists, for an equilibrium point, a positive definite function of type Fk) we say that 

the equilibrium point is stable with respect to the set of conserved quantities {Fx, Fk}. In 1958 G.K. 
Pozharitsky, see [5], had proved that it is sufficient to study the function ||(-Fi, Fk)\ \ in order to decide 
if an equilibrium point is stable with respect to the set of conserved quantities {Fx, ...,Fk}, (see |llj). 
Another method to decide if an equilibrium point is stable with respect to the set of conserved quantities 
{Fx, ...,Fk}, see [llj . is given by an algebraic method which reduces to study if the equilibrium point 
x e is isolated in the set of all the solutions of the algebraic system Fx(x) = Fx(x e ), ...,Fk(x) = Ff~{x e ). 
We also show that if the equilibrium point x e is not isolated in the set of all the solutions of the 
algebraic system given above, then it is impossible to construct a Lyapunov function in x e using the set 
of conserved quantities {Fx, ...,Fk}. 

We apply this algebraic method to decide the stability of an equilibrium point with respect to a set 
of conserved quantities for the problem of torque- free gyrostat. In Section 2 are presented some notions 
and results on Stability Theory which we apply in the study of our example. 

In Section 3 we present the mathematical model of a torque-free gyrostat and we give a set of two 
functionally independent conserved quantities. We find the set of uniform rotations and we first study 
their stability with respect to a single conserved quantity. In the cases when the vector of gyrostatic 
moment is situated along a principal axis of inertia of the gyrostat, we study the stability of an uniform 



od 
o 



i 



rotation with respect to the set of the two conserved quantities. We prove that an uniform rotation 
is stable with respect to the given set of conserved quantities if and only if it is stable in the sense 
of Lyapunov. It is interesting to see that exist some singular cases for which we cannot decide the 
Lyapunov stability of an uniform rotation using the algebraic method or the linearization method. 

In a future study we will apply the algebraic method which is used in this paper for the problem of 
the rotational motion of a gyrostat in the presence of an axisymmetric force field. We take advantage 
of the set of conserved quantities found in [3J. 

2 Lyapunov's direct method solving algebraic equations 

We consider an open set D C M™ and the locally Lipschitz function /:!)—>• M™ which generates the 
differential equation 

x = f(x). (2.1) 

We denote by x(-,Xq) the maximal solution of the above differential equation which verify the initial 
condition x(0,x o ) = x . A point x e S D is an equilibrium point of (|2.1[) if and only if f(x e ) = 0. An 
equilibrium point x e £ D is stable (or stable in the sense of Lyapunov) if for all e > there exists 
S > such that for all y in the ball B(x e ,S) and t > we have ||a;(t, y) — x e \\ < e (see [5]). The most 
important result for proving stability of an equilibrium point is given by Lyapunov's direct method. 

Theorem 2.1. Suppose there exists a continuous function V : D — > R satisfying the conditions: 

i) V(x e ) = 0; 

ii) V(x) > for x in a neighborhood of x e and x ^ x e ; 

Hi) t — > V(x(t,y)) is a decreasing function for all y G D. 

Then the equilibrium point x e is stable. 

A continuous function which satisfies the conditions i) and ii) is called a positive definite function 
in the equilibrium point x e . A continuous function V satisfying the hypotheses of the above theorem is 
called Lyapunov function in the equilibrium point x e . We introduce the following notion of stability. 

Definition 2.1. The equilibrium point x e of (12.11) is stable with respect to the set of conserved quan- 
tities {Fi, ...,Fk} if there exists a continuous function $ : R fc — > K such that x — > $>(F\, ....,Fj~)(x) — 
^(i 7 !, Fk)(x e ) is a positive definite function in x e . 

In the conditions of the above definition, the function x — > Fk)(x) — $(-Fi, Fk)(x e ) is a 

Lyapunov function in the equilibrium point x e . We have the obvious consequences. 

Theorem 2.2. Let x e be an equilibrium point and {Fx, F/~} be a set of conserved quantities for (12.11) . 

(i) If x e is stable with respect to the set {Fx, ...,Fk\ then it is stable in the sense of Lyapunov. 

(ii) Let q e {1, k} be an integer number. If x e is stable with respect to {F±, F q }, then it is stable 
with respect to {iq, Fk}. 

We have the following equivalent conditions for the stability of an equilibrium point with respect to 
a set of conserved quantities. 

Theorem 2.3. Let x e be an equilibrium point of (|2.ip and {Fx,...,Fk} a set of conserved quantities. 
The following statements are equivalent: 

(i) x e is stable with respect to the set of conserved quantities {Fi, F^}; 

(ii) x — > \ \(Fi, Fk)(x) — (F±, Fk)(x e )\\ is a positive definite function in x e ; 

(Hi) the system F\{x) = F\(x e ), ...,Fk(x) — Ff.(x e ) has no root besides x e in some neighborhood of x e . 
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In 1958, G.K. Pozharitsky had proved the equivalence between (i) and (ii), see [5], [TT] pp. 130. 
Equivalence between (ii) and (Hi) appears in [TT] pp. 151. In the paper [JJ, Aeyels had presented an 
interesting proof for the implication " (Hi) =>■ x e is Lyapunov stable" . 

The Theorem l2. 31 fm) gives an algebraic method for establishing Lyapunov stability of an equilibrium 
point. Moreover, it also shows that if the equilibrium point x e is not isolated in the set of solutions for 
the algebraic system of equations then it is impossible to construct a Lyapunov function in x e using the 
set of conserved quantities {F%, ...,Fk}. We will apply this algebraic method to study the stability of 
uniform rotations for a torque- free gyrostat. 

Using the implicit function theorem we have the following necessary but not sufficient condition for 
positive definiteness of the function given in Theorem 12.31 (ii), see |llj pp.151. 

Theorem 2.4. Let x e be an equilibrium point of (|2.1I) and ...,Fk} be a set of C 1 conserved quanti- 
ties. A necessary condition for the stability of x e with respect to set of conserved quantities {F%, Fk} 
is that the jacobian matrix ^LEl^iDsl ( Xe ) be of rank strictly smaller then k. 

For the case of one conserved quantity, i.e. k = 1, we have the well known result. 

Theorem 2.5. Let x e be an equilibrium point of (|2 . 1[) and F a conserved quantity. The following 
statements are equivalent: 

(i) x e is stable with respect to the conserved quantity F; 

(ii) x e is a strict local extremum of F. 



3 The stability of the uniform rotations of a torque-free gyro- 
stat 

For the problem of torque-free gyrostat we find the set of uniform rotations and we study their stability 
with respect to a conserved quantity. In the cases when the vector of the gyrostatic moment is situated 
along a principal axis of inertia of the gyrostat, we study the stability of an uniform rotation with respect 
to the set formed by two conserved quantities. Except two singular cases the Lyapunov stability problem 
for the free-torque gyrostat can be decided using the algebraic method with two conserved quantities 
and the linearization method. In the singular cases we decide the Lyapunov stability by studying the 
dynamics in an invariant set. 

The equation for the rotation of a torque- free gyrostat is given by (see [3"].|12]) 

id = (luj + ft) x C3, (3.1) 

where uj is the angular velocity, and I is the inertia tensor and fl is the constant vector of gyrostatic 
moment. We denote by 1%, I2 and ^3 the principal moments of inertia and suppose that 1% > I2 > I3. If 
we use the angular momentum vector M — luj then the equation becomes 

M = (M + ft) x \- l M. (3.2) 
It is easy to see that for the above dynamic we have two conserved quantities 

F 1 = ~M ■ l- l M, F 2 = ~(M + fl) ■ (M + ft). 

Next, we find the set of the uniform rotations. In the paper |10) was considered the differential equation 

N = N xI^N + axN, (3.3) 

where ael 3 . This equation is equivalent with the torque- free gyrostat equation (|3.2[) where a = — 
and making the change of variable M = N — ft. According to [10] the equilibrium points of (|3.3[) are of 
the following types: 

i. Nt = (0,0,0); 
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«• ^ = (A> A- A) fOT A G R\{£, h T 3 h 

iii. 7V 3 = (a, ^) if /xx = and ael; 

iv. iV 4 = (^ii2-,a,^_) if M2 = Oand a eM; 

v. N 5 = (j±h.,j^j- 2 ,a) if Ms = and a e K. 

Consequently, the uniform rotations of the torque-free equation (13. 2j) arc of the types: 
i- Mi = (-//!, -/i 2 , -/x 3 ); 

a- ^ = (l^Mi, T^M2, t^Ms) for A 6 R\{£, i £}; 

iii. M 3 = (/3, 7^M2, T^ft) if Mi = and /3 e M; 

iv. A? 4 = (j^f-iuP, t^jiIH) if M2 = and /3 e M; 

v. Ms = ( 1 £ K pi, 1 ±j;it i ,P) ifM3 = 0and/?GM. 

Analogous considerations are made in |12j . pp. 78-80, for finding the uniform rotations of the system 



First we study the stability of an uniform rotation with respect to one conserved quantity. 
Theorem 3.1. For the uniform rotations of a torque-free gyrostat we have: 

(i) The unique uniform rotation which is stable with respect to Fx is (0,0,0). This uniform rotation 
is of type M% obtained for A = 0. 

(ii) The unique uniform rotation which is stable with respect to F2 is M\ — (— /ii, — M2, — ^3)- 

Proof, (i) The uniform rotation (0,0,0) is the unique strict local extremum of the conserved quantity 
Fx. Using the Theorem 12.51 we obtain the result. 

(ii) The uniform rotation (— /Xi, — —^3) is the unique strict local extremum of the conserved quantity 
i<2- Using the Theorem 12.51 we obtain the enounced result. □ 

The uniform rotations found in the above theorem are the only uniform rotations which Lyapunov 
stability can be proved by using only one of the conserved quantities. For the rest of the uniform 
rotations it is necessary to consider both conserved quantities. Next, we study the stability of the 
uniform rotations with respect to the set conserved quantities {Fx,F2}. By direct calculus we obtain 
that the set of uniform rotations coincide with the set where the iacobian matrix d ( Fl 'f 2 ^ has the rank 
strictly smaller 2 and consequently, the necessary condition of Theorem 12.41 is fulfilled. 

In what follows we restrict ourselves to the cases for which the vector of gyrostatic moment /2 is 
situated along a principal axis of inertia of the gyrostat. 

3.1 The case /i 2 = A*3 = 

In this case we have the following types of uniform rotations: 

i. M1-2 = (q, 0, 0) where q E R; 

ii. Mi = (5^/^1,5,0) where q e R*; 

iii. M 5 = (j^fxx,0, q) where qeR*. 

First, we study the solutions of the algebraic system 

Fx(M) = Fx(M e ), F 2 (M) = F 2 (M e ), 
where M e is an uniform rotation. The above system of algebraic equations has the form: 

( M%_ Ml _ M t 2 e M| e M| e _ M£ 

< 12 + h ~ ii + 12 + is h (3.4) 
1 Mf + Ml = M? e + Mf e + Mf e - Ml - 2/uMi 
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where the unknowns are Mi,M 2 and M 3 . The system has at least the solution (Mi e) M% e , M3 e ). We 
want to see if this solution is isolated in the set of all the solutions of the algebraic system. For our 
study is preferable to change the variable Mi with x = M\ — M\ e . The algebraic system (|3.4p becomes 

f Ml A/f _ Af| e M| e x 2 +2xM 1<s 

) i 2 i 3 ~ i 2 "T" i 3 h (3.5) 

[ M| + M| = Mf e + M| e - (x 2 + 2xM le ) - 2m (x + M le ) 

with the unknowns x,M2 and M3. The system has at least the solution (0, M2 e ,M3 e ). The solution 
(Mi e , M2 e , M3e) of (|3.4j) is isolated in the set of all the solutions of this system if and only if the solution 
(0, M2e, Mfe) of (|3.5[) is isolated in the set of corresponding solutions. If we use the unknowns M| and 
Mf , then we have a linear system. Using Cramer's rule we can find the solutions of (|3.5[) . 

I. The uniform rotation of type Afi_ 2 . The solutions of (I3.5P verifies 

I M| = - / 2 ) + - h) + 1 J 



LI. If q = — j , then the system ([376]) becomes 

A ,f2 _ x 2 I 2 (Ii-I 3 ) 

M 3 ~ h(I 2 -I a ) \2 X i^ - 7 2j + 7,-/3 Ml J 

By our hypotheses we have I\ > 1% > I3 and if (x, M 2 ,M 3 ) is a solution of (|3.7p . then M| < 0. We 
deduce that (0,0,0) is the unique solution of the above system and consequently, it is isolated in the 
set of all the solutions. 

1.2. If q — — 1^\ 2 , then the system Q3.6[) becomes 



M 2 _ x^hjh-h) 
M 3 - Ji(J 2 -/ 3 ) 



(3.8) 



For |x| sufficiently small we have 



/ 2xl 2 (1 h(h-h) \\ f . . . 

sgn [- Il(l2 _ h) ~ h) + h _ h Mi ) ) = -sgn(x) • sgn(^i). 

For every \x\ sufficiently small such that sgn(x) = — sgn/xi we obtain a solution of (|3.8p and consequently, 
we have that (0, 0, 0) is not an isolated solution in the set of all the solutions. 

1.3. The case when q ^ — f^} 2 and q =/= — j^\ 3 ■ For x sufficiently small the terms in the righthand 
side of the system Q3.6P have the properties 



sgn Q 35 ^ 1 ~ ^ + q ( Tl ~ ^ + = s S n ( x ) ■ sgn{q{h - h) + h^i)- 

If sgn(g(/i — J3) + ■ sgn(q(I\ — I2) + hv>i) > 0, then exists r > such that a solution of the form 

(x, M2, M3) which verify x ^ has the property |x| > r. In this case, the solution (0, 0, 0) of the system 
p.6p is an isolated solution in the set of all the solutions. 

If sgn(q(7i — I3) + IxHi) ■ sgn(q(Ii — I2) + I1H1) < 0, then for every \x\ sufficiently small we have 
that the solutions of the system (|3 .6[) are of the form (x, M%, M3). We obtain that the solution (0, 0, 0) 
of the system l|3.6p is not isolated in the set of all the solutions. 



II. The uniform rotation of type M4. In this case the system (|3.5p is equivalent with the following 

system 

"'2 = T - / 1 (7 2 -7 3 ) x '" t- TT~r : 

j2 

'3 _ Ii(h-h) 
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As before, for every |x| sufficiently small we have a solution of the above system which is of the form 
(x, M 2 ,M 3 ) and consequently, (0, q, 0) is not isolated in the set of all the solutions of (|3.9[) . 

III. The uniform rotation of type M5. The system (|3.5|) is equivalent with the following system 

( m 2 - _£!MW3) 

The solutions of the above system are (0, 0, q) and (0, 0, —q). The solution (0, 0, q) is isolated in the set 
of all the solutions of (|3.10p . 

Summarizing, we obtain the following result. 

Theorem 3.2. For a vector of gyrostatic moment along the first axis of inertia we have the following 
stability results. 



(i) An uniform rotation of type M\_ 2 = (q, 0, 0) is stable with respect to the set of conserved quantities 
{F 1 ,F 2 } if and only if q G (-oo,-^) U [-^,00) and M i > or q e (-00, U 

r J1U1 



j^fj- , 00) and /ii < 0. 

(ii) An uniform rotation of type M4, = ( j _j Ml; 9j 0) with q ^ is not stable with respect to the set of 
conserved quantities {F\ 1 F 2 \. 

(Hi) An uniform rotation of type M5 = ( j^j Hi, 0, q) with q ^ is stable with respect to the set of 
conserved quantities {Fi,F 2 }. 

We prove that the uniform rotations which are not stable with respect to the set {Fi , F 2 } are 
Lyapunov unstable. In the paper [10] is proved that the equilibrium points of the system (13.3[) have the 
properties: 

(i) for A G (j-, j-) an equilibrium point of type N 2 is spectrally unstable; 

(ii) an equilibrium point of type -ZV4 is spectrally unstable when a =/= 0. 
Consequently, we have: 

(i) for q G (-77^77, - j^j^ ) and fix > or for q G {— f^) 3 , — j^f^ ) and ^1 < an uniform rotation 
of type Mi_2 is spectrally unstable and consequently, it is unstable in the sense of Lyapunov; 

(ii) an uniform rotation of type M4 is spectrally unstable and also it is unstable in the sense of 
Lyapunov. 

The Lyapunov stability or instability of M e = {— ,0,0) cannot be decided using the set of 
conserved quantities {Fi,F 2 } or using the linearization method. This uniform rotation is spectrally 
stable and it is not stable with respect to the set of of conserved quantities {Fi,F 2 }. The instability in 
the sense of Lyapunov of this uniform rotation will be proved by studying the dynamics on the invariant 
set 

M = {M\ Fi(M) = Fi(Me), F 2 (M) = F 2 (M e )}. 

Theorem 3.3. The uniform rotation (— , 0, 0) is unstable in the sense of Lyapunov. 

Proof. The projection of the vectorial differential equation (|3.2j) on the first axis using the variables 
x, M 2 and M3 is 

x = (^ - hM 2 M 3 . 

By using (|3.8p we have 

2x\l x -h) (\ h{h-h) 

X = TTFt FT9 ( o^v 1 — 3j H r r Ml 



r((h-h) 2 \2 yi J/ h-h 

First we consider the case [i\ > 0. Suppose that we have x(0) > — jj^j^frj^r^j Hi and M 2 (0)Mz(0) < 
0. Consequently, we obtain that x(0) < 0. In this case there exists t* > such that x(t*) = 
— (ii-i2){ii-i 3 ) Mi w hich implies that our uniform rotation is unstable. In the case fix < we have 
analogous considerations. □ 
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Remark 3.1. In this case an uniform rotation is stable with respect to the set of conserved quantities 
{Fx, F 2 } if and only if it is stable in the sense of Lyapunov. 

3.2 The case fii = fx 3 = 

We have the following types of uniform rotations: 

i. Mi_ 2 = (0,<Z,0) where q E R; 

ii. A% = (q, j^V2, 0) where q E R*; 

hi. M 5 = (0, j^fJ. 2 ,q) where q E R* . 

Using the method of the previous section and by analogous calculations we obtain the following result. 

Theorem 3.4. For a vector of gyro 'static moment along the second axis of inertia we have the following 
stability results. 

(i) An uniform rotation of type Mx- 2 = (0, q, 0) is stable with respect to the set of conserved quantities 
{F 1 ,F 2 } if and only if q G jf^] and li 2 >0 or q E [7^ , - andfi 2 <0. 

(ii) An uniform rotation of type M3 = (q, j^zj- ^ 2 , 0) with q ^ is stable with respect to the set of 
conserved quantities {Fx,F 2 }. 

(Hi) An uniform rotation of type M5 = (0, 7^377/^2,9) with q ^ is stable with respect to the set of 
conserved quantities {Fx,F 2 }. 

In the paper |10) is proved that an equilibrium point of type N 2 is spectrally unstable for A S (y-, j-). 

Consequently for q E R\[-j^, and fi 2 > or q E R\[j^, -7^77] and fi 2 < an uniform 

rotation of type Mi_2 = (0, q, 0) is spectrally unstable and also it is unstable in the sense of Lyapunov. 

In this case the stability (in the sense of Lyapunov) can be decided using the stability with respect 
to the set of conserved quantities {Fx, F 2 } and the linearization method. An uniform rotation is stable 
with respect to the set of conserved quantities {Fx , F 2 } if and only if it is stable in the sense of Lyapunov. 

3.3 The case fix = fi 2 = 

We have the following types of uniform rotations: 

i. A?i_ 2 = (0, 0, q) where qER; 

ii. M 3 = (q, 0, fi 3 ) where qER*; 

hi. M 4 = (0, q, j^r^) where qER*. 

Theorem 3.5. For a vector of gyro static moment along the third axis of inertia we have the following 
stability results. 

(i) An uniform rotation of type Mx-i — (0, 0, q) is stable with respect to the set of conserved quantities 
{Fx,F 2 } if and only if q E (-00, U 00) and ^ > or q E (-00, j^) U 00) 
and us < 0. 

(ii) An uniform rotation of type M3 = (g,0, 7^77/^3) with q ^ is stable with respect to the set of 
conserved quantities {Fx,F 2 }. 

(Hi) An uniform rotation of type M4 = (0, q, 7^7-^3) with q ^ is not stable with respect to the set of 
conserved quantities {Fx,F 2 }. 

In the paper [10] is proved that the equilibrium points of the system (|3 - 3[) has the properties: 

(i) for A E (j-, j-) an equilibrium point of type N 2 is spectrally unstable; 

(ii) an equilibrium point of type N4 is spectrally unstable when a^0. 
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These results implies: 

(i) for q E ( j*yf ) j^t 3 ) an d > or for g <E ( f ^| , 7^^) and /i 3 < an uniform rotation of type 
Mi_2 is spectrally unstable and consequently, it is unstable in the sense of Lyapunov; 

(ii) an uniform rotation of type M4 is spectrally unstable and also it is unstable in the sense of 
Lyapunov. 

The Lyapunov stability or instability of M e = (0, 0, jfzj^) cannot be decided using the set of con- 
served quantities {F\, F 2 } or using the linearization method. This uniform rotation is spectrally stable 
and it is not stable with respect to the set of of conserved quantities {Fi,F 2 }. The instability in the 
sense of Lyapunov of this uniform rotation will be proved by studying the dynamics on the invariant set 

M = {M\ Fx(M) = Fi(M e ), F 2 (M) = F 2 (M e )}. 

The proof is analogous to the proof of the Theorem 13.31 

Theorem 3.6. The uniform rotation (0,0, j^zj^) * s unstable in the sense of Lyapunov. 

An uniform rotation is stable with respect to the set of conserved quantities {F\, F 2 } if and only if 
it is stable in the sense of Lyapunov. 
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